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Abstract 

At very high energies, pair production (7 — > e~^e~) exhibits many interesting features. The 
momentum transfer from the target is very smah, so the reaction probes the macroscopic properties 
of the target, rather than individual nuclei. Interference between interactions with different atoms 
reduces the pair production cross section considerably below the Bethe-Heitler values. At very 
high energies, photonuclear interactions may outnumber pair production. 

In contrast, in crystals, the interaction amplitudes may add coherently, greatly increasing the 
cross sections. Pair production in matter-free magnetic fields is also possible. The highest energy 
pair production occurs at high energy particle colliders. This article will compare pair production 
in these very different regimes. 



1 



INTRODUCTION 



Pair production was first observed by Anderson in 1932, who used the process to discover 
the positron Q . Our modern theoretical understanding of pair production stems from work 
by Bethe and Heitler Q . At high photon energies, the pair production cross section reaches 
a constant value 0, 

where Xq is the target radiation length, in units of column density (mass/length^, usually 
g/cm^), A is the atomic mass of the target (usually in g/mole), and A^^ = 6.022 x 10^^ is 
Avogadro's number ^. Xq is also commonly given in units of length; this length can be 
multiplied by the density p to give Xq in column density. 

Eq. n applies at high photon energies, > 10 GeV, when the momentum transfer from 
the target to the pair is sufficiently small. The atomic electrons screen the nucleus from the 
incident photon/pair, limiting the effective impact parameter to the Thomas- Fermi radius. 
In this complete screening limit, the cross section is independent of photon energy. The 
differential cross section is usually given in terms of x = E/k, where E is the electron 
energy j^: 
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This lowest order cross section is symmetric in x and 1 — a;, the electron and positron energy. 

At high energies, two regimes are of interest. Most conversions produce pairs near thresh- 
old, with pair mass Mp ^ 2m, where m is the electron mass. The momentum transfer 
required to produce a pair decreases linearly as the photon energy rises. The longitudinal 
momentum transfer required to produce a lepton pair with energies E and k — E is 

,„=t-v/(I3B)^-Vl^«^^^«f (3) 

We take c = 1 throughout this review. The approximation in the last two formulae cover 
are accurate for the vast majority of the cross section, but fail for very asymmetric pairs, 
where x — > or x 1. 

As k rises and q\\ drops, the uncertainty principle requires that the location of the conver- 
sion becomes more and more delocalized. The pair production is delocalized over a region 
known as the formation length: 
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The subscript '0' shows that this is for pair conversion in free space; future sections will 
consider modifications to // in dense media. The formation length grows linearly with the 
photon energy. When the formation length is larger than the typical inter-nuclear separation, 
interactions with different atoms are no longer independent. 

This loss of independence leads to many interesting consequences. In amorphous ma- 
terials, the pair production cross section decreases as the photon energy rises, invalidating 
Eqs. H and |21 The conversion probability is no longer proportional to the number of target 
atoms. In moderately thick targets, the conversion probability may even be proportional 
to the logarithm of the thickness! The reduced cross sections appear in a number of 'real 
world' situations. For example, they affect electromagnetic showers produced by very high 
energy cosmic rays. 

In crystals, the long formation length allows pair conversion amplitudes from regularly 
spaced atoms to add coherently. The strong magnetic fields can greatly increase the pair 
production cross section. This enhancement may be useful in accelerator design. 

For very-high- mass pairs, the momentum transfer from the pair is high, and production of 
these pairs can probe short distance scales. High-mass pairs have been used to test quantum 
electrodynamics (QED). Neglecting the transverse momentum of the pair (a relatively small 
effect), the pair mass is 

Ml = -f^. (5) 

Figure H shows the pair mass spectrum, da/dMp. The cross section is strongly peaked near 
threshold, x ~ 0.5 and Mp ^ 2m, In this regime, the electron and positron have a small 
(non-relativistic) relative velocity, and the cross section is enhanced. At even slightly higher 
masses, the differential velocity is large and the cross section scales as 

da 1 
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This l/M^ depend ence continues as long as the momentum transfer to the target is small. 
When the momentum transfer rises, screening is no longer complete, and the structure of 
the target becomes significant. In this regime, the cross section decreases faster than 1/Mp. 

More recently, QED has been studied at particle colliders, which can reach shorter dis- 
tance scales. Pair production reactions such as e~^e~ — > e~^e~e~^e~ are an important channel 
for these studies. Another important channel for testing QED has been pair production 
at heavy ion colliders. Here, the very strong fields are of interest. For gold or lead, the 
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FIG. 1: da/dMp for pair production from a 100 GeV photon. The cross section is strongly peaked 
near threshold. Away from the peak, the cross section scales as 1/Mp. The spectrum is almost 
independent of the photon energy. For high enough pair masses, the interaction is sensitive to the 
targets internal structure, and the cross section drops falls off more rapidly. 

perturbative expansion coupling Za ~ 0.6, and it is not unreasonable to expect significant 
non-perturbative effects in heavy ion collisions. Here, a = /he ~ 1/137 is the electromag- 
netic coupling constant, with e the charge on the electron. 

High energy electrodynamics in matter has been previously reviewed by Akhiezer and 
Shulga ^, Baier and Katkov ||7|, and by Klein |5||. Although it does not address pair 
production, Ter Mikaelians book has a very interesting, largely classical discussion of 
high ener gy bremsstrahlung. Interactions in crystals have been considered by Palazzi 
S0rensen jlfl| . and by Baier, Katkov and Strakhovenko Pair production in heavy ion 



collisions has been recently reviewed by Gerhard Baur and his collaborators [1^ 

Section 2 of this review considers pair production from isolated atoms. Section 3 will 
discuss pair production in bulk media. Section 4 will cover production in discrete fields and 
crystals. Section 5 will consider pair production at colliding beam accelerators, and Sec. 6 
will draw some conclusions. 
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FIG. 2: A schematic diagram of (a) direct pair production on a target with charge Z and (b) pair 
production via the Compton process. Diagram (a) is usually dominant, and the Compton process 
(b) is mainly visible via it's interference with direct production. 

PAIR PRODUCTION FROM ISOLATED ATOMS 

Pair production from an isolated atom occurs schematically via a two-photon process, 
as is shown in Fig. |2ta); In the target frame of reference, the incident photon fluctuates 
to a virtual e~^e~ pair. The pair combines with a virtual photon from the target to form a 
real e~^e~ pair. In the competing Compton process. Fig. I2tb), the photon is first absorbed 
by a nucleus or electron, which then emits a virtual photon, which decays to a e~^e~ pair. 
The Compton amplitude is relatively small, and the process is usually observed via it's 
interference with the direct reaction P]. This review will focus on the kinematics of pair 
production, which determines many characteristics of the reaction. 

As long as // > a, where the Thomas- Fermi atomic radius a = 0.8Z~^^^aQ (oq = 5.3 
nm is the Bohr radius) then the nascent pair couples to the target atom as a whole. This 
is the complete screening limit; at distances larger than a, the electrons screen the nuclear 
charge. In this limit, the cross section for pair production is independent of k. For heavy ions 
(heavier than iron), the limit is reached within 1% for k > 10 GeV; for lighter ions, higher 
photon energies are required. For hydrogen, the finite-energy correction is still slightly larger 
than 1% for k = 100 GeV 

For light ions, the radiation length Xq scales as 1/Z^ - the interaction amplitude depends 
on Z, and the cross section on Z^. However, as Z rises, higher order diagrams with multiple 
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photons may be important. The most important (largest) higher order terms account for 
the fact that the pair is produced in the potential well of the target nucleus; the other terms, 
due to fluctuations in the incoming photon, are small. Bethe and Maximon calculated the 
corrections due to the potential by solving exactly the Dirac equation for an electron in a 



Coulomb field 



14| . The solution used Furry wave functions in parabolic coordinates. These 



radiative corrections reduced the cross section by a constant amount 



li 



where rg = 2.8 fm is the classical electron radius and 

For a heavy atom like lead, /(82) = 0.33, and the Coulomb correction reduces the cross 
section about 10%, essentially independent of k. 

For the differential cross section, the leading order correction depends only slightly on x. 
However the next term is large for a; ^ or a; ^ 1 . The sign of the second term differs 
for the two extremes in x, introducing an asymmetry between the electron and positron 
momentum spectra. These asymmetries have also been studied with perturbation theory; 



for heavy nuclei , sig nificant asymmetries have been predicted theoretically 
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and observed 



experimentally |l8|. 

Tsai made detailed calculations of the screening factors and radiative corrections for 
different materials. He also considered pair conversion on the atomic electrons in the target. 
The complexities of these calculations can be easily hidden by adjusting Xq. With this, for 
Z > 4 Q, 

Xo= (^l^^[z2[ln(184Z-i/3)_/(z)] + Zln(1194Z-2/3)]^ (9) 

The first logarithmic term is for nuclear scattering, while the second is for scattering from 
atomic electrons; f{Z) is the Bethe-Maximon correction. For light nuclei with Z < 4, the 
logarithmic terms are inaccurate, and Z— dependent constants are usually used. With this 
Xo, Eqs. HI and El remain accurate and easily usable. At 10 GeV, these cross sections agree 
with the data to within 2% [iQ^ . 

At very high energies, other channels may contribute to the cross section. For example, 
the virtual photon in Fig. |21can be replaced with a virtual boson. Although the amplitude 



for the mediated reaction is not large, the interference with the purely electromagnetic 
process can introduce ~ 5% asymmetries between the electron and the positron, especially 
for large Mp 

There is a paucity of data on pair production from photons with A; > 10 GeV. The high- 
energy frontier has moved to e'^e" and hadron colliders, as will be discussed in Section 4. 
There is been considerable data on pair production in crystals, but the effects of the crystal 
overshadow any subtleties in the individual photon-atom interactions; crystals are discussed 
in Section „ 



PAIR PRODUCTION IN MEDIA 



As the photon energy rises, the formation length becomes large enough that it encom- 
passes more than one atom. Then, the amplitudes for interacting with multiple atoms must 
be added, and interference between interactions with different atoms needs to be considered. 

Landau and Pomeranchuk considered this problem in 1953 
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They used classical 



electrodynamics to study bremsstrahlung. The emission of photons with momentum k and 



energy k = \k\, 
solid angle Q is 



rom an electron starting with velocity vi and ending with velocity V2 into 
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k ■ Vi — k k ■ V2 — k 

This equation is usually applied to an interaction of an electron with a single atom. Lan- 
dau and Pomeranchuk realized that, when atoms are close enough together, the interactions 
cannot be separated. The emission depends only on the initial and final velocities, and is 
independent of any intermediate state velocities. The scattering due to nearby atoms adds 
together. 

As long as the angle 6^ between vi and V2 is less than I/7, where j = E/m is the Lorentz 
boost of the electron, then the bremsstrahlung radiation is proportional to \vi — V2\^ and 
the radiation from the summed independent (with uncorrelated directions) scattering is the 
same as if the radiation from each scatterer was added independently. In this case, we define 
Av = Vi — V2, the electron momentum change is g = 7m | Av\ , and 



dkdVL 
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Integrating over the solid angle gives the total emission: 

dl 2Z^eY 



(12) 



dk Svrm^ 

Equation [T^ holds as long as the momentum transfer g < m. At larger momentum transfers, 
the radiation is reduced. The cross section for classical bremsstrahlung can be found from 
Eq. EH 

Landau and Pomeranchuk calculated the radiation from a long electron trajectory that 
was determined by multiple scattering, with the scattering spread evenly along the path. 
They determined the interference from radiation from different points on the path. Radiation 
from nearby points added in-phase, while at larger separations, the phase was random. The 
phase coherence holds for shorter and shorter distances as the multiple scattering increased. 
For energetic electrons radiating low energy photons, the reduction in coherence reduces the 
total emission. 

The conditions for reduced emission can be estimated from the multiple scattering in a 
formation length. When the multiple scattering angle within If is larger than I/7, then the 
multiple scattering reduces the coherence, and, with it, the radiation. In a distance //, the 
mean multiple scattering angle is ^ 

- §^ (13) 



where Es = ^-n ja = 21.2 MeV. An additional multiplicative term, 1 + 0.038 In (Zj/Xq) 
is sometimes included to account for the non-Gaussian nature of Coulomb scattering. This 
non-Gaussian nature is particularly apparent for thin targets, where the number of scatters 
is small. However, since Eq. is only used for back-of-the-envelope calculations, this 
multiplicative factor will not be used here. 

Neglecting the non-Gaussian scattering tails, and requiring 6ms < I/t? the emission of 
photons with energy k from electrons with energy E is suppressed when 

k < —TTT- (14 

In this limit, the bremsstrahlung spectrum is altered from the Bethe-Heitler prediction 
da /dk ^ 1/k to a harder spectrum, da/dk ^ l/\/k\ 

Although pair production is not a classical process. Landau and Pomeranchuk noted that 
similar reasoning should apply for it. By using crossing symmetry to relate bremsstrahlung 
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and pair production, they predicted that for photons with very high energies, the pair 
production cross section scales as 

This is very different from the Bethe-Heitler energy-independent cross section. At high 
enough energies, photons become penetrating particles. 

This cross section can be derived using the uncertainty principle Q] and relying on the 
proportionality between the pair cross section and formation length. The number of atoms 
in the formation zone is proportional to //. For coherent interactions, the probability scales 
as the number of atoms squared. The pair conversion amplitudes add in-phase over the 
formation length, so the interaction probability scales as /j, and the probability per unit 
length scales as Z/. Multiple scattering reduces //, and with it the cross section. 

Multiple scattering changes the electron and positron direction, but not their momenta, 
reducing their longitudinal velocity. The longitudinal velocity reduction Av is 

^ = [1 - oosiOMs)] - (16) 

The slowing reduces the longitunal momentum of the electron and positron. The excess 

momentum is transferred to the medium during the scattering. Including the multiple 

scattering from the positron and electron, 6ms+ and 0ms-, 

^ m'k Eels+ . jk - E)ehs- ^ rn'k E%k 
^" 2E{k-E) 2 2 2E{k - E) 2E{k - E)Xq ^ ' 

When 6ms > V7 ^^e electron or positron, the multiple scattering terms dominates. This 

happens when 

2 

Ttl 

ho > ^Xo, (18) 

i.e. when 

In this regime, the multiple scattering angles depend on If, which itself depends on 9ms- 
This leads to a quadratic equation for If. When the multiple scattering dominates. 



a _ If _ / kEiPM Mp Elpm 



O'er Ifo V — E) m V k 
Here, Elpm is a material dependent constant. 



(20) 



Elpm = -TE^f ~ 7.7 TeV/cm ■ Xq. (21) 



For lead, Elpm = 4.3 TeV. Sometimes Elpm is defined differently, usually differing by 
a numerical factor (often 2). Since Elpm depends only on Xq, the pair production cross 
section may be easily calculated for mixtures. 

Pair production is suppressed for photons with k > E^pm, i-e. in the TeV region and 
above. In the strong suppression limit {k ^ Elpm), the formation length scales as ~ v^, 
and the cross section is suppressed as /////o ~ l/Vk. Suppression is largest for low- mass 
pairs. This is not surprising; the larger Mp, the shorter I^q (Eq. H}, and the less of an 
opportunity for multiple scattering to contribute to q\\. 

While the pair production cross section falls at high energies, the photonuclear cross 
section rises. Photonuclear interactions occur when the photon fluctuates into a quark- 
antiquark {qq) pair, which then interacts hadronically with a target nucleus. Very high 
energy photons may also interact directly with quarks in the nucleus 125^ . 

When the photon energy is high enough, photonuclear interactions may dominate over 
pair production [2^ . Fig. El compares the pair production cross section in lead and water 
(using Migdals suppression, see Sec. 1^ with the photonuclear cross sections. The photonu- 
clear cross sections have not been measured at these energies, but should increase slowly 
with energy. 

For lead, the photonuclear cross section, ahad is about 15 mb when the photon-nucleon 
center of mass energy, y/s^, is less than 100 GeV, and 15(A/s^7v/100GeV)°'^ mb at higher 
energies. For water, ahad = 1-9 mb for ^/s^ < 100 GeV, and 1.9(a/ s^/v/100GeV)°'^ mb at 
higher energies. These parameterizations roughly follow Fig. 5 of Ref . (2^ ; the oxygen cross 
section is scaled from carbon as the nuclear cross section, A^^^, and the water cross section is 
the sum of its constituent cross sections. Any estimate in this energy range requires a large 
extrapolation, so there is considerable uncertainty. The energy-dependence in Ref. is 
reasonably conservative. 

There are also significant uncertainties in the pair production cross section. The cross 
section calculations may fail when the suppression becomes very large (order of a ~ 1/137). 
When LPM suppression is very large, higher-order reactions such as — » e'^e~'yA are more 
important. With their larger inherent q\\, these reactions are less affected by by multiple 
scattering. 

With these caveats, for both materials the photonuclear cross sections are larger than 
the pair production cross sections for k > A x 10^^ eV. The crossover energy is similar for 
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FIG. 3: Cross sections for 7 — > e~^e~ (solid line) compared with the cross section for photonuclear 
interactions (dotted line) in (a) lead and (b) water. Caveats for the curves are discussed in the 
text. 

both materials. For lighter materials, ahad/<^BH is larger, while for heavier materials, LPM 
suppression is larger. The two effects largely cancel out. 

Even with the uncertainties, the trend is clear: at very high energies, photonuclear in- 
teractions become more important, and may even dominate over pair production. Even 
at a few 10^^ eV, electromagnetic showers may develop significant hadronic components, 
complicating their identification. 

Production of fi^^^ and r+r^ pairs, 7 l^l^ is reduced by a factor [m/miY, where mi 
is the lepton mass. So, neither process is very important here. Of course, in some parts of 
phase space (such as for high mass pairs), the suppression is smaller. 

One other suppression mechanism should be mentioned: bremsstrahlung can also suppress 
pair production p, A produced leptons can emit a bremsstrahlung photon and lose 
momentum. When the bremsstrahlung occurs in the pair production formation zone, the 
combined reaction is a single, higher-order process, the bremsstrahlung will increase q\\, 
decreasing If, and with it, the pair production cross section. 

In a crude model, the kinematics changes when If > Xq. This only happens at very 



11 



high energies, when the LPM effect is already important. Including LPM suppression, 
If > Xq when k > E^Xo/h, or, for lead, k > 4 x 10^^ eV. This is only approximate; 
while bremsstrahlung suppresses pair production, the reverse is also true: pair production 
suppresses bremsstrahlung. The suppression of one reaction reduces its impact on the other 
reaction |2fi] and a combined calculation is needed to properly understand pair production 
and bremsstrahlung in this energy regime. 

The energies required to study LPM suppression of pair production are well beyond 
the reach of existing accelerators. However, LPM suppression is important to studies of 
very high energy cosmic rays (including neutrinos). When high energy cosmic rays hit the 
atmosphere, they interact, forming air showers which develop into billions of lower energy 
particles. Cosmic rays with energies up to 3 x 10^° eV have been observed jt28i|. If the incident 
particles are photons (unlikely, as discussed in Section then LPM suppression reduces the 
cross sections and thereby lengthens the shower, increasing the amount of energy that reaches 
the ground. Proton induced showers are much less affected, but still have an electromagnetic 
component that is subject to the LPM effect. However, the suppression is only relevant for 
the highest energy showers {29]. As will be discussed in Section „ cosmic ray showers have 
been used to study the LPM effect. 

Very high energy cosmic ray z/g, z/^ and Ur may be visible at next generation neutrino 
detectors, like IceCube 30]. Neutrino interactions produce high energy leptons along with 
hadronic showers from the struck nucleus, z/g showers are most affected by the LPM effect, 
because of the very high energy electron. However, electrons from r decay and the electro- 
magnetic component of the hadronic shower may also be affected by LPM suppression. 



Thin Targets and Surface Production 

In very thin targets, the formation length may be larger than the target thickness T. 
When this occurs, the target interacts as a single unit. Including LPM suppression. If > T 
when 

2mT 

k> ^2.5TeV -Tifim). (22) 

When If > T, LPM suppression may be less than in a very thick target. 

For very thin targets, the multiple scattering in the target is negligible, and the Bethe- 
Maximon cross section is retained, independent of k. The Bethe-Maximon cross section 
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applies when the q\\ from the multiple scattering (the second term in Eq. El with T replacing 
If) is less than that due to the pair production (the first term in Eq. [T7|) . This occurs for 
target thicknesses T < {m/ EsYXq. 

For slightly thicker targets, {m? / E'^)Xq < T < If, the conversion probability depends 
on the total scattering in the target. This regime has been studied only for bremsstrahlung 



j31[. Assuming that the behavior is the same, for k S> Elpm, the conversion probability 
scales as In {EgT/m^Xo) - as the logarithm of the target thickness! 

When T > If, the target no longer acts as a single unit. However, there is an increased 
probability of interactions within If of the target surface. The additional interactions are 
known as transition pair creation, analogous to transition radiation. The plasma frequency 
of the medium, Up plays a key role. For solids, hup is 30-80 eV, so about hoUp ^ 10~'^m. 
Taking a = hup/ m, j32l | 

27r^ 35 256 ^' ^ ^ 

independent of the incident photon energy. In a typical solid a ~ 10""^ and -P(a) ~ 10~^°, 

so this effect is generally negligible. 

At higher energies, there are additional surface effects due to multiple scattering. A 

new type of transition radiation occurs as the electrons electromagnetic fields rearrange 

themselves to account for the multiple scattering in the medium. Ternovskii calculated 

the differential probabilities for bremsstrahlung and. pair production in plates of moderate 



thicknesses 



33[. For plates that are at least Tq = aXo/27r^, (1 < ^ < 2 rises slowly with k 



to account for the increasing shower length), the differential probability for pair production 



IS 

dN a 



The scaling with target thickness T is far from linear! This equation only holds for relatively 
thin targets. For thick targets, surface transitions are treated as an add-on to production 
in the bulk of the target. Ternovskii gave formulae for surface transitions which add-on to 
Migdals calculation of bulk transitions discussed in the next section. 



Migdal Formulation for Pair Production 

In 1956, Migdal published the first accurate calculation of suppression of pair production 
{3^ . He focused on bremsstrahlung, but also provided formulae for pair production using 
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crossing symmetry. He represented the transition probabilities with a density matrix, and 
then averaged the density matrix over the possible positions of the scattering centers. This 
led to a Boltzman transport equation for the radiation. He solved this with the Fokker- 
Planck technique. 

This approach used the quantum mechanical matrix elements, with proper account of 
electron spin, photon polarization and nuclear recoil. The scattering was realistically (ran- 
domly) distributed in space (albeit with a Gaussian distribution for angles). 

Migdal expressed his results in terms of two functions G{s) and 0(s), where s depends 
on the photon and electron energies, and on the total radiation (Migdal used s to discuss 
bremsstrahlung, and s for pair production; the two variables are closely related): 



^8Eik-E)m' ^ ' 

The factor ^(s) increases from 1 to 2 as the pair production moves from no suppression 
(s > 1, the Bethe-Heitler limit) to strong suppression (s <^ 1): 

m = 1 > 1) 

e(s) = l + ln(s)/ln(si) (si < s< 1) 

m = 2 {s< s,) (26) 

with Si = Migdal did not discuss mixtures of materials. However, the radiation 

is dominated by high— Z materials, it is not too far wrong to take the Z of the heaviest 
atom in a mixture. Alternately, one could, following the approach used for Xq take a 
weighted average, weighting the different materials by their abundance times Z^. 

The equations for s and ^ can be solved recursively. However, ^ (s) varies very slowly, and 



Stanev and collaborators found an approach that avoids the recursion 
following Eq. with ^ = 1. Then, 



3- 



They defined s' 



as') = 1 (^' > 1) 

e(s') =2 {s'< V2s,) (27) 

where /i = In (s')/ In {\/2si). 
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FIG. 4: Migdals G{s) (dashed line) and (t>{s) (solid line). 



The differential cross section for pair production is 



da 



(28) 



dE SXoNAkK"''^' ' L k2 

Migdal gave infinite series for 0(s) and G{s). However, they may also be represented as 
polynomials 



(f){s) = 1 — exp 



6s[l + (3 - 7r)s] + sV(0.623 + 0.79s + 0.658s^) 



(29) 



and 



■^{s) = 1 - exp[-4s - 8sV(l + 3.96s + 4.97^2 - O.OSs^ + 7.5s^)] (30) 

with G(s) = 3i/j{s) — 2(f)(s). Figure IH shows these functions. 

The limit s — oo, ^(s) = 1 corresponds to low photon energies, where 0(s) = Gil) = 1 
and the Bethe-Heitler cross section holds. 

The limit s — 0, ^(s) =2 corresponds to strong suppression. Then, 0(s) = 6s, G(s) = 
127rs^ and the cross section is 



da 



4A 



0(s). 



(31) 



dE 3XoNAk V 

Figure |31 shows the pair production cross sections vs. E/k in lead for different photon 
energies. At very high energies, most conversions produce asymmetric pairs. 



15 




0.25 -\ 



I I I " I 

0.25 0.5 0.75 1 

x=E/k 

FIG. 5: Differential cross section for pair production, Eq. 1281 in a lead target, for different photon 
energies. Cross sections are shown for photons with k = 1 TeV (top), 10 TeV, 100 TeV, 1 PeV, 
10 PeV, 100 PeV and 1 EeV (bottom). Here, n is the number of atoms per unit volume. When 
E{k — E) 1.2kELPM, the product (j){s)^{s) rises slightly above 1 and the conversion probability 
for high energy photons is slightly larger than at smaller k. These curves can be used for other 
materials by scaling k by E lp m i^ead) / Elpm', the bottom curve also applies for 3.2 EeV photons 
in iron (E^pm = 13.6 TeV). 

Figure El gives the suppression (compared to Bethe-Heitler) as a function of pair mass 
in lead for photons with energies from 44 TeV to 440 PeV. The suppression is largest for 
small Mp. The Bethe-Heitler da/dMp shape is independent of k, so the cross section is 
the product of this suppression factor with the cross section from Fig. ^ The bulk of the 
cross section is near threshold, so the overall suppression is dominated by the measurement 
around 2m ^ 1.12 MeV. This calculation neglects the lepton transverse momenta, but most 
of the consequent inaccuracy should cancel for the suppression. Very roughly, in the strong 
suppression limit, for Mp ^ 2m the cross section scales as da/dMp ^ ^/M^. 

When the cross section suppression is large, the penetration depth of electromagnetic 
showers increases significantly. Not only is the total cross section reduced, but either the 
electron or positron takes most of the energy. 
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FIG. 6: The suppression factor S = ctlpm / <^bh as a function of pair mass for photons with energies 
44 TeV (top curve), 440 TeV, 4.4 PeV, 44 PeV and 440 PeV (bottom curve). The smah dips below 
S = 1 around Sm?k/ Eipm are side-ripples to the dips in Fig|SJ These curves can be applied 

to other materials by scaling k by ELPM(}ead)/ELPM] the bottom curve also applies for 1.36 EeV 
photons in iron {Elpm = 13.6 TeV). 

Newer Calculations 

Several recent papers have presented more sophisticated calculations of LPM suppression. 
These calculations consider targets as an integral whole, and therefore naturally handle finite 
target thicknesses. However, many authors considered only bremsstrahlung. Often, too little 
information is given for an outsider to easily apply the calculations to pair production. They 
also do not consider mixtures of atoms with different Z. Blankenbecler and Drell used an 
eikonal approach to study bremsstrahlung in finite targets js^. 

Zakharov used a light cone path integral approach to study bremsstrahlung. He used a 
transverse Green's function to represent e^e~ pair scattering from the atom and used it to 
solve the Schroedingers equation. He used an imaginary potential for the Greens function, 
proportional to the cross section for an e^e~ pair to scatter from the atom. 

In the strong suppression limit, the scatterer is a dipole with length p; the scattering 
cross section cr(p) = C(p)p^ where C(p) varies slowly with p js?!. With a constant C(p), 
this is Gaussian scattering. For an infinitely thick target, the results are almost identical 
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to Migdals. With a slowly varying C{p) to model Coulomb scattering, the Zakharov cross 
section has more energy dependence than Migdals ^{s)- At higher energies, Zakharov should 
find slightly lower cross sections than Migdal. 

For regimes with moderate suppression, Zakharov adopted a more accurate Coulomb 



potential than a dipole, with separate potentials for the nucleus and it's electrons. 



?his 

3- 



separation is important for low-Z atoms where the electrons play a significant role 
With these details, Zakharov found a good fit to the SLAC E-146 data, except for the 
carbon. 

V. N. Baier and V. M. Katkov made one of the few explicit calculations of pair pro- 
duction The framework was similar to Migdal, but with a quasiclassical eikonal 
approximation for the scattering. This quasiclassical approach has been criticized as being 
inapplicable in the regime where it was used because the electron impact parameter cannot 



be treated classically 



41[. However, Baier and Katkov argue that, when the angular mo- 



. The eikonal approach allowed them to 



mentum is large, their approach is acceptable 
include Coulomb corrections in the potential, an important factor in improving the accuracy 
of the calculations. They included a separate term in the potential to account for scattering 
from the atomic electrons. In the no-suppression limit, the authors obtain the fully-screened 
Bethe-Maximon Coulomb-corrected cross sections. 

Baier and Katkov note that, for heavy elements, the Coulomb correction is larger than 
the uncertainties in the SLAC E-146 data [7]. The E-146 experimental analysis normalized 
the Migdal calculations (without radiative corrections) to a radiation length that included 
Coulomb corrections, effectively multiplying the two effects together. This was the limit 
of the then-current technology. Although the results matched the data, this may not have 
been an optimal approach. 

When the suppression is strong, Baier and Katkov give a polynomial expression for the 
total pair production cross section 



ctbm V k 



1 _ i.672W%^ - 2.192^ + 7^(1^ + 0.274) 

\l k k 4L/ AElpm ' 



(32) 



where Li = 183Z~^/^e~-^, with / from Eq. |H1 Here, asM is the Bethe-Maximon cross section 
(apparently mis-labelled as Bethe-Heitler in Eq. 2.35 of Ref. jQ). The polynomial in 
Eq. 1221 reduces the cross section at moderate photon energies. However, at large enough 
energies, the polynomial terms drop out, but the last correction term in Eq. |221 increases 
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the cross section to 2-10% above the Migdal result. 



Experimental Measurements of Suppression 

Suppression of pair production can only be observed for photon energies k > Elpm- 
Since Elpm >1 TeV, even for the densest materials, it has not been possible to study this 
suppression at accelerators. Some studies have used cosmic rays, but the statistics are limited 
by the high energy cosmic ray flux. The earliest studies considered low-energy electrons in 
photon induced showers. Fowler, Perkins and Pinau studied 47 electromagnetic showers 
with an energy above 1 TeV in an emulsion stack, measuring the distance between the 
initial conversion point and the first daughter pair downstream of the primary |^|. At high 
photon energies, the average distance rose, as expected for decreasing cross sections. The 
statistics were limited because few sufficiently energetic photons penetrate the atmosphere. 
To improve the statistics, some later experiments flew emulsion detectors in high-altitude 
balloons where the flux of high-energy photons is higher. They achieved somewhat 
better statistics, but with qualitatively similar results. 



One of the most precise studies was by Kasahara in 1985 j45|. He studied the development 
of electromagnetic showers with energies of order 100 TeV, in lead-emulsion chambers. At 
these energies, LPM suppression significantly increases the penetration depth of the showers. 
After removing the contamination from hadronic showers, Kasahara found that the shower 
profiles were consistent with Migdals cross sections, but not the Bethe-Heitler formulae. 

Because of the limited statistics possible in cosmic ray studies, the best studies of the LPM 
effect have used accelerators to study the suppression of bremsstrahlung. Bremsstrahlung 
can be studied with electron beams with energies E <^ Elpm, as long as photons with 
energies k < E^/E^pm are studied. 

The first accelerator experiment used 40 GeV electrons generated from 70 GeV protons at 
the Serpukhov U-70 accelerator Bremsstrahlung photons from carbon, aluminum, lead 
and tungsten targets were detected in a sodium-iodide calorimeter. The electrons were bent 
away from the calorimeter with a magnet. The group studied photons with energies between 
20 and 80 MeV. At higher photon energies, the LPM effect was unimportant, while at lower 
energies, the backgrounds were very high. Even in the signal region, there were significant 
backgrounds from bremsstrahlung in air and scintillation counters, and muon contamination 
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in the beam. By taking ratios of their data (lead/aluminum and tungsten/carbon), they 
were able to estimate the degree of suppression. They found suppression larger than Migdal 
predicted, albeit with large errors. 

A later experiment at the Stanford Linear Accelerator Center, SLAG Experiment E-146, 
made precise measurements of bremsstrahlung from 7 targets, from carbon to uranium, 

3 HQ- 



in beams of 8 and 25 GeV electrons 



between 200 keV and 500 MeV. The photons were detected in a position-sensitive segmented 
BGO calorimeter 50 meters away from the target. The small solid angle and position 
sensitivity helped reduce backgrounds, especially from synchrotron radiation. Data was 
taken at 2 different settings of calorimeter gain (i.e. photomultiplier tube high voltage) 
to maximize the dynamic range. The electron beam-line was kept entirely in vacuum to 
eliminate bremsstrahlung from air or target windows. 

The group used a tertiary electron beam, but, because the primary beam was also elec- 
trons, beam contamination was small. The beam was usually run at an average of 1 electron 
per pulse, 120 pulses/second; events containing zero or multiple electrons were rejected 
during analysis. 

The calorimeter was calibrated using cosmic rays and a 500 MeV electron beam. The 
group also used the downstream magnet for an electron spectrometer, measuring the electron 
energy loss. Because of the limited resolution for electron energy loss, this was useful only 

ClS db check. 

Figure U\ shows the E-146 data for the aluminum targets. Targets with thicknesses of 
3% and 6% of Xq (3.12 mm and 5.3 mm) were used. Events containing a single electron 
were selected, and the photon energy spectrum measured. The photon energy is plotted in 
logarithmic width bins, with 25 bins per decade of energy. Then, the approximate Bethe- 
Heitler photon spectrum dN/dk ^ 1/k transforms to dN/d\n{k)/XQ ^ 0.13 - a straight 
line. The probability for a single electron to undergo two independent interactions in the 
target tilts the spectra slightly. This was dealt with by using a Monte Garlo simulation. 

The data shows significant suppression. For k < 2 x 10~^i?, an additional suppression 
mechanism is required to explain the data. This second mechanism involves the produced 
photon, which also interacts with the medium, through forward Gompton scattering. Glas- 
sically, this is described via a dielectric constant e 7^ 1. 

The group has also studied thin targets, where T < If for relevant photon energies. 
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For thin targets, the radiation was above the Migdal calculation, and consistent with the 
calculations that treated the entire target as a single radiator. 

More recently, an Aarhus-CERN-Florence-Free University, Amsterdam-Cape Town col- 
laboration studied bremsstrahlung from 149-287 GeV electrons from iridium, tungsten and 
tantalum targets 5^ . They used a CERN test beam. Figure |H1 shows their tantalum target 
data plotted in a manner similar to the E-146 data. With the higher electron energies, 
they observed suppression up to k/E 0.16, a significant fraction of the spectrum. The 
goal of this experiment was to directly observe the overall increase in radiation length as 
the radiation decreases. They observed roughly a 30% decrease in radiation at the highest 
electron energy, consistent with expectations from Migdals calculation. 



The collaboration fitted their spectra to measure E^pm for their targets [5l|. For the 
densest target, iridium, their measured Elpm averaged 1.97 ± 0.16 TeV, in quite good 
agreement with the theoretical prediction of 2.25 TeV. However, for the lighter targets, 
tantalum and copper, the measured Elpm were considerably below the expected values - by 
40% for copper. The reasons for the difference are not clear, but the collaboration suggests 
that possibly bremsstrahlung from atomic electrons could be unsuppressed. The experiment 
also uses data from a carbon target for a normalization spectrum; if the carbon spectra were 
somehow compromised, this data could affect their measured Elpm- It may be worth noting 
that the E-146 carbon-target data is poorly described by calculations based on Migdals work. 

PAIR PRODUCTION IN CONSTANT FIELDS AND CRYSTALS 

Over 50 years ago, H. Robl noted that for very high energy photons, matter is not nec- 
essary for pair production: even static electromagnetic fields suffice j^. This is analogous 
to synchrotron radiation. The photon conversion depends on the parameter 

-1^ 



where Bj^ is the magnetic field perpendicular to the photon direction and = iv?/eh = 
4.4 X 10^^ Gauss is the critical field. The probability of photon survival for a distance d is 

P{d) = exp (-ad) (34) 

where js^ 

<.= ^A-3(2/3x) (35) 
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and Ki/3{x) is a modified Bessel function and Ae = 386 fm is the Compton wavelength of 
the electron. This equation is for the constant field approximation, and may fail for rapidly 
varying fields. Figure El shows the attenuation coefficient a for photons of different energies. 
The absorption is small for x ^ 1, or A; < 21713^/ Bj_. In this limit a = 0.46 exp (— 4/3x)- 
The absorption rises sharply as a function of k for x > 1, reaching a peak for k = 12mi?c/ B±. 
At this maximum, a ^ 0.1{a/Xe)Bj_/Bc. As x continues to rise, the attenuation coefficient 
drops slowly; for x > 1000, a ^ 0.6%^^^^ ~ B'^/^k~^^^. The attentuation increases with 
increasing B±, but, for fixed B±, drops slowly as k rises. 

There are some minor discrepancies between various studies. Erber 0] uses a dimension- 
less auxiliary function T(T) to calculate a. His Eq. 3.3d, defines T(T), with numerically 
matching small and large T limits in Eq. 3.3c. However, Fig. 9 and Table VI of Ref. 53| 
give T(T) are half as large as the equations; the equations appear to be correct. Ref. 5J| 
plots photon attenuation following the equations in Erbers review. 

The momentum transfer to the field required for pair production is the same as for an 
atomic target, Eq 01 Additional terms are required to account for the magnetic bending 
which, like multiple scattering, reduces the electron and positron longitudinal velocity. The 
field bends electrons smoothly, with 

Ob = 5%^. (36) 

in contrast to the 9 ^ found for multiple scattering. Including the magnetic bending 
for the electron {9b-) and positron {9b+), the longitudinal momentum transfer is 

Ml E9l_ , {k-E)9l^ Ml e^Bll)k 
^11" 2A; + 2 + 2 ~ 2k ^ 2E{k - E) ^'^^^ 

This leads to a a quartic equation for //. At low photon energies, magnetic bending is small 
and If rises linearly with k. It reaches a maximum, and then decreases . For % ^ 1, 

( hE{k-E) \ _ ( hkm^ \ 

so ^ k^l'^/B'^f - the same scaling found from the more detailed calculations. 

Since B^. is so large, this process is only important for photons with quite high energies. 
In fact, conversion in an external magnetic field has not been observed. However, it does 
have important astrophysical applications. One example involves very high energy cosmic 
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rays. Very high energy photons may pair convert in the earth's magnetic field 
For a k = 10^° eV photon perpendicular to the earth's field, the attenuation length is about 
100 km. The earth's magnetic field extends to altitudes of lOOO's of kilometers, far above the 
atmosphere where lower energy showers originate. The produced electrons interact with the 
earth's magnetic field and emit synchrotron radiation. The combination of pair production in 
the field and synchrotron radiation leads to an electromagnetic shower, even in the absence 
of matter. The field-induced shower largely stops developing when the average particle 
energy drops so that x ^ 1 (~ 10^^-^ eV). The produced particles propagate downward 
until they reach the atmosphere, where the shower develops conventionally. The initial 
high-altitude reactions gives these showers a head-start. This early development is helpful 
in differentiating photon initiate showers from proton initiated showers. Only a handful of 
showers with e nerg ies above 10^° eV have been observed, but the primary particles cannot 
all be photons j28|. 

Pair conversion can also occur in the strong magnetic fields that surround many astro- 
physical objects. For example, many neutron stars support surface magnetic fields around 
10^^ G. In these fields pair conversion occurs quickly for A; > 50 MeV. 

Although pair conversion from an external magnetic field has not been seen in the labo- 
ratory, a fairly close analog has been studied: pair conversion in crystals. When a photon 
traverses a crystal parallel to one of the axes, the magnetic fields of the individual atoms 
add together to act as a single, extremely strong field. The effective fields are largest when 
the photon direction is well aligned with a crystal axis, within a characteristic angle given 
by the atomic lattice spacing divided by the formation length. For 100 GeV photons, this 
angle is a few mrad. Often, the target is mounted in a goniometer which allowed precise 
adjustments of the crystal orientation. Pair conversion is studied as the crystal orientation 
is varied; when the axis is properly aligned, pair production is enhanced j60|. 

Two factors affect the enhancement. Impurities and other lattice defects lead to uneven 
atom spacings which reduce the coherence; high quality crystals are required. The quality 
requirement limits the choice of target material. At the highest photon energies, variations 
in atomic spacing due to lattice vibrations limit the coherence. Sometimes, crystals are 
cooled to reduce these vibrations. The second factor is the divergence (angular size) of the 
photon beam. In a beam with a large divergence, many photons will not be within the 
characteristic angle of the crystal axis. 
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The simplest manifestation of the strong coherent fields is channelling; charged particles 
in crystals may be strongly deflected. Channeling studies show that the effective fields can 
reach 25 Megagauss in silicon crystals js^]. Crystal channeling can in some cases supplant 
conventional accelerator magnets, particularly for beam extraction js^ . 

With the increased bending, production of bremsstrahlung photons is also greatly 
enhanced. This technique is mature enough that plans are advanced to use coherent 
bremsstrahlung in a diamond crystal to generate high-energy polarized photon beams for a 
series of experiments js^. 

The fields can also cause photon conversions. At CERN, pair conversion in germanium, 
tungsten and iridium crystals has been studied with photons with energies up to 150 GeV. 
The pair production cross section was enhanced over the isolated atom Bethe-Heitler predic- 
tion; the enhancement grew linearly with energy, reaching a factor of 7 for k = 150 GeV IgI]. 
The total and differential (with respect to E/k) cross sections agreed with the predictions 
for coherent enhancement. 

When a magnetic field also contains matter, both fields must be considered together. The 
effects of the nuclei and the fields must be combined to find the conversion probability. This 
is a complex calculation. Here we examine one limit, the reduction in pair conversion on a 
target due to the motion caused by a moderate (x < 1) magnetic field. Here, the magnetic 
field plays a similar role to multiple scattering. When the magnetic bending terms in Eq. 
l37l dominate, the magnetic field suppresses the bremsstrahlung. The suppression equation 
is quartic, and, in the limit of strong suppression, the cross section is 

f kEB 

where 

Eb = ^ (40) 

plays a similar role to E^pm- This approach only applies for very high energy photons 
in relatively low-density matter. At high densities, LPM suppression dominates. When 
X > 1, pair conversion in the magnetic field is important. A parallel analysis finds a similar 
suppression for bremsstrahlung. The bremsstrahlung calculation is supported by detailed 
calculation using kinetic equations {g^ . 
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PAIR PRODUCTION AT PARTICLE COLLIDERS 



Pair production also occurs at particle colliders. The colliding particles each emit a 
photon; the two photons collide to form an e^e~ or other lepton pair. Lepton pair production 
at e^e~ colliders has been studied intensely over the past 30 years. These studies have been 
used for a variety of tests of QED {g^. The reactions have also been studied as backgrounds 
to other interesting processes, especially those that are beyond the standard model. They 
are also used to measure the internal (electromagnetic) structure of the photon |6J]. 

Although these processes follow the same diagrams as pair production, Fig. 1, there are 
some significant differences from the fixed target regime. The reaction e^e~ e~^e~e~^e~, 
studied at e^e~ colliders, involves many more diagrams than appear in Fig. 1. Second, 
the photons are virtual. For e~^e~ colliders, and ep colliders, the virtuality (effective mass) 
is large, while in proton and ion colliders it is smaller. This review will focus on photo- 
production with almost-real photons (those with small virtuality) at ep, pp/pp and heavy 
ion colliders. Finally, the symmetric initial state leads to some interesting effects. There 
is a vast literature on photoproduction at colliders; here we only consider a few selected 
examples to illustrate the similarities and differences with fixed target reactions. 

Muon and electron pair production, ep epe^e~ and ep ep^^fi^ were studied by the 
HI collaboration at the HERA ep collider jci^. They observed lepton pairs with masses up 
to 80 GeV, where electroweak corrections to QED are significant. The study was done to 
search for physics beyond the standard model, but the collaboration found good agreement 
with the standard model. 

In colliders, the ions are fully stripped of electrons. There is no screening, and, as 
the beam energy rises, photoproduction occurs at larger and larger impact parameters. 
Consequently, the cross section continues to rise with energy indefinitely. For example, in 
gold on gold collisions at a center of mass energy of 200 GeV per nucleon at Brookhaven's 
Relativistic Heavy Ion Collider, the cross section is expected to be 33,000 barns, rising 
to 200,000 barns in 5.5 TeV per nucleon lead on lead collisions at CERN's Large Hadron 
Collider Pair production can occur at ion- ion separations up to several fim. 

Because Za ~ 0.6, higher order corrections to the cross section could be large, and 
pair production may probe non-perturbative strong-field QED. This possibility encouraged 
many theoretical efforts. One approach extended the Bethe-Maximon approach to collisions 
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of heavy ions. The effect of the Coulomb potentials of both ions was to reduce the cross 
section by 25% at RHIC and 14% at the LHC 3- 

Tony Baltz made an all-orders calculation of pair production in ultra-relativistic heavy 
ion collisions. He found an exact solution to the time-dependent Dirac equation jo^. In an 
appropriate gauge, the Coulomb potential from an ultra-relativistic particle simplifies into 
a simple two-dimensional potential, and pair production can be calculated exactly. Surpris- 
ingly, his results matched the lowest order perturbative result ^67^ , without any Coulomb cor- 
rections. However, a later, more careful calculation found results that included the Coulomb 
co„.ect,o„s 0. The theoretical ,ss„e. a.e discussed e.ewhere i. this vo.un,e, so they w,U 
not be considered further here. 

In most existing data on e^e~ production in ion collisions, at least one of the ions is light 
(typically sulfur). Experiment NA-45 at the CERN SPS studied 200 GeV per nucleon bare 
sulfur ions hitting a platinum target and found cross sections consistent with the lowest 
order perturbative result 0|- The measured cross section in a restricted kinematic region 
is 0.99^*^0 28 of that expected from lowest order QED. Here, the errors are combined in 
quadrature. This result could be used to put limits on some models of radiative corrections. 
Other SPS studies found similar results, with general agreement to lowest order calculations 



3, 



III. 



At ion colliders, the charges are large enough that multiple interactions between a single 
ion pair are very possible. For example, two pairs may be created by a single ion-ion reaction 
1721 1. Multiple pairs have not been studied yet. However, the STAR collaboration has studied 

_ nn 

e^e pair production accompanied by mutual Coulomb dissociation of the two nuclei |73l. l74| . 
The mutual Coulomb dissociation criteria selects collisions where the impact parameter 
is less than about 30 fm, considerably smaller than for unselected pair production. The 
reduction in impact parameter leads to stronger colliding fields, and therefore to increased 
higher order corrections. Despite this, the collaboration found that their results were well 
described by lowest order QED. They did find that it was necessary to include the photon 
virtuality in the calculation to explain the pt spectrum of the pair. 

Multiple interactions are most common at small impact parameters, where the field den- 
sities are strongest. For this reason, they are a good place to search for higher order effects. 
Despite this, the STAR data are consistent with lowest order QED. Because of the restricted 
impact parameter range, this data cannot be directly compared with existing studies of ra- 
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diative corrections 

The absence of atomic electrons allows some new reactions, such as pair production with 
capture, where an electron is produced bound to one of the participating nuclei. This reaction 
has a number of notable applications, most notably the first production of antihydrogen. 
Antihydrogen is produced when positrons were produced via pair production, bound to 

n 

antiprotons 75]. At heavy ion colliders, the cross section for electron capture is sizable. 



of the order of 100-200 barns at RHIC and the LHC [7^. The single-electron atoms are 
lost from the colliding beams, reducing the luminosity, and giving pair production some 
importance for accelerator design. 



CONCLUSIONS 



Pair production at high energies covers a varied palette of physics. High mass pairs have 
been a sensitive test of QED, and pair production at e~^e~ and ep colliders continues to be a 
laboratory to search for new physics. These colliders allow the highest energy tests of QED. 

Lower mass pairs have long formation lengths. These reactions are distributed over 
many atoms, and probe bulk characteristics of the material. When the photon energy is 
large enough, pair production is suppressed below the Bethe-Heitler cross section, and many 
'obvious' scaling laws fail. At very high photon energies, the pair production cross section 
may be smaller than the photonuclear cross section, and photonic showers will look like 
hadronic showers. 

Pairs may also be produced in strong magnetic fields. This has significant implications 
in astrophysics. These strong magnetic fields also appear in crystals when the photon is 
aligned with one of the crystal axes. These aligned photons may coherently convert to e^e~ 
pairs. 
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discussions. Kai Hencken, Gerhard Baur, Roman Lee and Hristofor Vankov all provided 
useful suggestions. This work was supported by the U.Sy. Department of Energy under 
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FIG. 7: SLAG Experiment E-146 data (points) compared with Monte Garlo predictions using 
the Bethe-Heitler (dashed Hne), LPM suppression (dot-dashed hne) and with LPM and dielectric 



suppression (soUd Hne). The simulations include transition radiation. Adapted from Ref. 



49|. 
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FIG. 8: Bremsstrahlung spectrum from (a) 287 GeV, (b) 207 GeV and (c) 149 GeV electrons 
striking a 0.128 mm thick iridium target. The photon energy scale is logarithmic, with 25 bins per 
decade if energy. The dotted line is a simulation using the Bethe-Heitler spectrum, while the solid 



line includes LPM suppression. From Ref. 



5C|. 
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FIG. 9: Attenuation coefficient (inverse attenuation length) for photons in a magnetic field with 
perpendicular components of 0.35 Gauss (the Earth's average field - solid line), 100 Gauss (short 
dashes), 10^ Gauss (long dashes) and 10^ Gauss (dotted line). 
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